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XXVIII. Dynamical Problems regarding Elastic Spheroidal Shells and Spheroids of 
Incompressible Liquid. By Professor W. Thomson, LL.D., F*B.S. 

Eeceived August 22, — Bead November 27, 1862. 



1. The theory of elastic solids in equilibrium presents the following general problem. :— 
A solid of any shape being given, and displacements being arbitrarily produced or 
forces arbitrarily applied over its whole bounding surface, it is required to find the dis- 
placement of every point of its substance. The chief object of the present communica- 
tion is to show the solution of this problem for the case of a shell consisting of isotropic 
elastic material, and bounded by two concentric spherical surfaces, with the natural 
restriction that the whole alteration of figure is very small. 

2. Let the centre of the spherical surfaces be taken as origin, and let at, y, z be the 
rectangular coordinates of any particle of the solid, in its undisturbed position, and #+#? 
y+$, 2+7 the coordinates of the same particle when the whole is in equilibrium under 
the given superficial disturbing action. Then, by the known equations of equilibrium 
of elastic solids, we have 
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m—j^n and n denoting the two coefficients of elasticity, which may be called respectively 
the elasticity of volume, and the rigidity. A demonstration of these equations, with defi- 
nitions of the coefficients, will be found in & 71 of an Appendix to the present commu- 
nication. 



3. For brevity let 
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so that 5 shall denote the cubic dilatation at the point (#, y, z) of the solid. Also, for 

d z d 2 d 2 
brevity, let the operation ^+^o+^ be denoted by y 2 . Then the preceding equa- 
tions become 



n^y+mj^O. 

j 
4k2 



(3) 
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4. In certain cases, especially the ideal one of an incompressible elastic solid, the 
following notation is more convenient : — 

p the mean normal pressure per unit of area on all sides of any small portion of the 
solid, round the point #, y, z. Then (below, § 21) 

i>=-(—M(£+f+S); w 

and the equations of equilibrium become 



m—^n dx 

%V 2 /3 * t=0, 

m—fn dy 
m-—\ndz 



• ».■»• 



• (5) 



5. If the solid were incompressible, we should have m=oo and 

«?# dfi dy n 
dx^ay^dz ' 

which must be taken instead of (4), and, along with (5), would constitute the four diffe- 
rential equations required for the four unknown functions on, j3, y, p *. 

d d 

6. To solve the general equations (3) or (5), take -r of the first, j~ of the second, 



and j- of the third ; and add. We have thus 



(%+m)V 2 ^=0, (6) 

or, which is in general sufficient, 

V 2 &=0 (7) 

If, now, an appropriate solution of this equation for & is found, the three equations (3) 
may be solved by known methods, the first of them for a, the second for /3, and the third 
for y, — the arbitrary part of the solution in each case being merely a solution of the 
equation v 2 w=0. These arbitrary parts must be determined so as to fulfil equation (2) 
and the prescribed surface conditions. 

The complete particular determination of $ cannot, however, in most cases be effected 
without regard to a, )3, y ; and the order of procedure which has been indicated is only 
convenient for determining the proper forms for general solutions of the equations. 

7. First, then, to solve the equation in $ generally, we may use a theorem belonging 
to the foundation of Laplace's remarkable analysis of the attraction of spheroids, which 
may be enunciated as follows. 

If the equation V 2 &=0 is satisfied for every point between two concentric spheres of 

* See Professor Stokes's paper " On the Friction of Fluids in Motion, and the Equilibrium and Motion of 
Elastic Solids," Cambridge Philosophical Society's Transactions, April, 1845. 
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radii a (greater) and a! (less), the value of I for any point of this space, at distance r from 
the centre, may be expressed by the double series 

V +V X +V 2 + &c. 

+ v>- i +v>- 3 +v;r- 5 + &c, 

of which the first part converges at least as rapidly as the geometrical progression 

/ \ 2 / \ 3 

r /r\ I r 
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and the second at least as rapidly as 

V v 

r 

— if V|, Vi denote homogeneous functions of %, y, z of the order % each satisfying, con- 
tinuously, for all values of #, y 9 z, the equation 

V 2 V=0. 

A proof of this proposition is given in Thomson and Tait's 'Natural Philosophy,' chap. i. 
Appendix B. It is also there shown, what I'believe has been hitherto overlooked, that 
V„ V-, as above defined, cannot but be rational and integral, if i is any positive integer. 

8. To avoid circumlocution, we shall call any homogeneous function of (#, y, z) which 
satisfies the equation V72V— ft 

a "spherical harmonic function," or, more shortly, a "spherical harmonic.'' Thus V 4 
and V-, as defined in § 7, are spherical harmonics of degree or order i; and V-r"" 21 *"* 1 , 
being also a solution of V 2 V=0, is a spherical harmonic of degree — (i+1). We shall 
sometimes call the latter a spherical harmonic of inverse order i. Thus ii { being any 
spherical harmonic of integral degree % and therefore necessarily a rational integral 
function of this degree, u^ 21 " 1 is a spherical harmonic of degree — (i+1), or of inverse 
order L 

If we put— (£+l)==y\ and denote this last function by fy, then we have 

and thus it appears that the relation between a spherical harmonic of positive degree i 
and of negative degree j is reciprocal. The general (well known) proposition on which 
this depends is that if V 4 is any homogeneous function of (#, y, z) of degree i 9 positive 
or negative, integral or fractional, Y^- 2 *- 1 is also a solution of the equation V 2 V=0 
(see Thomson and Tait's 4 Natural Philosophy,' chap. i. Appendix B.). 

A spherical harmonic of integral whether positive or negative degree, satisfying the 
differential equation continuously for all values of the variables, will be called an " entire 
spherical harmonic," because such functions are suited for the solution of acoustical and 
other physical problems regarding entire spheres or entire spherical shells. 

A spherical harmonic function of (x, y, z) will be called a " spherical surface-har- 
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monic " when the point (#, #, z) lies anywhere on a spherical surface having its centre 
at the origin of coordinates. A spherical surface-harmonic is therefore a function of 
two variables, angular coordinates of a point on a spherical surface. If Y 8 - denote such 
a function of order i, positive and integral, then Y/* and Y^"*" 1 are what we now call 
simply spherical harmonics; but sometimes we shall call them, by way of distinction, 
"spherical solid harmonics." Functions Y,., or spherical surface-harmonics of integral 
orders, have been generally called " Laplace's coefficients " by English writers. 

9. From the theorem enunciated in § 7, we see that the general solution of our 
problem, so far as $ is concerned, is this : — 

*=%zT(Vi+Yr*- 1 ) (8) 

10. Now because the equation V 2 u=0 is linear, it follows that differential coeffi- 
cients of any solution, with reference to #, y> 2, or linear functions of such differential 
coefficients, are also solutions, Hence the terms V,. and V^" -2 '- 1 , of &, give harmonics of 

the degrees i— 1 and — (^+2), in — -, — , ~~. To solve equations (3) we have therefore 
only to solve 

where <p n denotes an entire spherical harmonic of any positive or negative degree, n, 
Trying 

which is obviously the right form, we have 



V 2 m=A< 



h*+<-£-»*+'£)*+«™} 



But, because <p n is a homogeneous function of #, #, z of degree w 3 



(*i+y i+4y= n *»-> 



and because it is a spherical harmonic, 

V 2 <p M =0. 
We have also 

V 2 (r 2 )=6, 
by differentiation. Hence 

V 2 %=A.2(2w+3)0„, 

and therefore the complete solution of the equation 

V 2 u=<p n 

n 2(2n + 3) rn ' 

where V denotes any solution of the equation 

V 2 V=0. 
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11. Hence, by taking for <p n the terms of -j- , T , T above referred to (§ 10), and 

(tX u/U QiZ 

giving n its proper value, e— 1, or --(£+2), for each term as the case may be, we find, 
for the complete solution of (3), the following: — 
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where ^, «4 v { , v^ w^ w'i denote six harmonics, each of degree i. 

12. But in order that these formulae may express the solution of the original equations 
(1), the functions u, v, &c. must be related to the functions V so as to satisfy equations 
(2) and (3). Now, taking account of the following formula, 
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V dx J* dy\ dy J 
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which becomes simply 



2ty 4 , 



if fy is a spherical harmonic of any degree i (whether positive or negative, integral or 
fractional), we derive from (9) by differentiation, and selection of terms of order % and 
of order inverse i (or degree — ^— 1), 



dx*dy*dz~ ^ 

where, for brevity, we put 
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Hence, to satisfy (2) and (8), 
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13, Using these in (9), we conclude 
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J 



for a complete solution of the general equations (1), the equations of equilibrium of an 
isotropic elastic solid. The circumstances for which this solution is appropriate will 
be understood when the general proposition of § 7 is duly considered. 

14. It remains to show how the harmonics u u v { , w { , u' i9 v' i9 w\ are to be deter- 
mined so as to satisfy the superficial conditions. Let us first suppose these to be 
that the displacement of every point of the bounding surface is given arbitrarily. 
Let 2Af , 2B f , 2Q be the harmonic series f , expressing the three components of the 
displacement at any point of the outer surface of the shell, and 2AJ, 2BJ, 2C». the 
corresponding expressions for the given condition of the inner surface. Thus the 
surface-equations of condition to be fulfilled are 



r=a < 



a: 

(3: 

L 7- 
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2 A,,) 
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4^i 5 J 



(13) 



where a and a 1 denote the radii of the outer and inner surfaces respectively, and 
Ai , B^ , C f , A[ , B| , Ci spherical surface-harmonics of the order L 

15. Now collecting from the series (12) of § 13, which constitute the general expres- 
sions for a, j3 5 7, those terms which, being either solid spherical harmonics of degrees 
% and — 4—1 , or such functions multiplied by r 2 , give, at the boundary, surface-harmonics 
of the order % and equating the terms of this order on the two sides of equations (13), 
we have 

* For the case i=0, the terms u._ v v._ v iv'.^ may be omitted; but their Ml interpretation would be to 

"" A 

express a displacement without deformation. Thus il_ v being of degree — 1, cannot but be — , where A is 

a constant ; and therefore u. r~ 2i+l becomes A when i=zO* 

f That is, series of terms each of which Is a spherical surface-harmonic of integral order i. That any 
function ; arbitrarily given over an entire spherical surface, may be so expressed, is a well-known theorem. 
A demonstration of it is given in Thomson and Tait's < Natural Philosophy/ chap, i. Appendix B, § s. 
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16. These six equations would suffice to determine the six harmonics %, v i9 w { , «4 V& wj, 
if $ i+l and ^_ x were known. For, since each of those six functions is a homogeneous 
function of #, y, z of order % each of them divided by r l is a function of angular coordi- 
nates relative to the centre, and independent of r ; and therefore if, for instance, we denote 
u . by r% and u'i by rW, we have two unknown quantities & and m r to be determined by 
the two equations of condition relative to cc for the outer and the inner surface. These 

equations may be written as follows, if we further denote -~±l by r'B> and ^ J^ 



by r"" 1 '"" 1 ^, because these are homogeneous functions of the orders £ and —v 
ively : 
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Resolving these equations for w and «r', and returning to the original notation instead 
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where, for brevity, 
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Introducing, also for brevity, the following notation, 
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we have the expressions for u t and u 4 given below. Dealing with the equations of 
condition relative to j3 and y, and introducing an abbreviated notation B i9 B^, C i9 CJ 
corresponding to (17), we find similar expressions for v i9 v i9 Wi , w i9 as follows: — 
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17. It only remains to determine the functions ^ and ^', which we can do by com- 
bining these last equations with (10) of § 12. Thus, changing i into i+l in (17) and 
into i— 1 in (18), applying equations (10) of § 12, and taking advantage of the following 
properties, 

v 2 ^. +2 =o, V 2 (^-)=0, &a, 
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These equations, used to determine the two unknown functions \|/ 4 and $, give 
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18. The functions $ and %J/ being expressed in terms of the data of the problem by 
equations (22), (23), (17), (18), (16), we have only to use (19) and (20) in (12) to find 
the following expression of the complete solution : — 
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19. This solution leads immediately, through an extreme case of its application, to 
the solution of the general problem for a plate of elastic substance between two infinite 
parallel planes: — Given the displacement of every point of its surface, required the 
displacement of any interior point. For if we give infinite values to a and a!, and keep 
#■— a! finite, the spherical shell becomes an infinite plane plate. 

20. It is, however, less easy to deduce the result in this way from the solution for 
the spherical shell, than to apply directly the general method of § 6 to the case of the 
infinite plane plate. We shall return to this subject (§ 31, below), when the details of 
the investigation will be sufficiently indicated. 

21. A very important part of the general problem proposed in § 1 remains to be 
considered, — that in which not the displacement, but the arbitrarily applied force, is 
given all over the surface. To express the surface-equations of condition for such data, 
we must use the formulae expressing the stress (or force of elasticity) in any part of an 
elastic solid in terms of the strain (or deformation) of the substance. These are 
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where P, Q, R are the normal tractions (which when negative are pressures) on the 
faces of a unit cube respectively perpendicular to the lines of reference OX, OY, OZ, 
and S, T, U the tangential forces along the faces respectively parallel, and in the direc- 
tions in these planes respectively perpendicular, to OX, OY, OZ (see Appendix, § 70). 
22. In terms of these we have the following expressions for the components F, G, H 
of the force on a unit area perpendicular to any line whose direction cosines are /, </, h : — 

F =If+TJg+Th, ' 

G=Uf+Qg+Sh, > . (26) 

lx=Tjr -j-S^ -f-BA 
4l2 
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(see " Elements of a Mathematical Theory of Elasticity," Philosophical Transactions for 
1856, p. 481). 

23. Using the expressions (25) in (26), we find 

^(— )(£+!+£)/+.(£/+$*£»)+•(£*■ &+£*). • <«> 

and symmetrical expressions for G and H. 

24. If now we suppose f 9 g, h to denote the direction-cosines of the normal at any 
point is i y, z of the surface of an elastic solid, the surface condition, when force, not 
displacement, is given, will be expressed by equating F, G, H respectively to three 
functions of the coordinates of a point in the surface, quite arbitrary except iii so far 
as they must balance one another in order that equilibrium in the body may be possible ; 
and therefore they must fulfil the following integral equations :— • 

jjFd[Q=0, jjGdO=0, $M2=0, (28) 

tf(Hy-<fe)^^ . , . (29) 

where dQ, denotes an element of the surface at the point (#, y, z\ and the double 
integrals include the whole surface of application of the forces F, G, H. 

25. For our case of the spherical shell, with origin of coordinates at its centre, we have 

f=7> ff=h h =7> W 

and the last triple term in the expression (27) for F may be conveniently written thus : — 



7 5 ® m * \ J 



Then, for brevity, putting 



and 



ax+Py+yz=Z, • ( 32 ) 



4+4+*£ =r l; (33) 

d 

where ~ prefixed to any function of #, y> z will denote its rate of variation per unit of 

length in the radial direction ; and using (2) of § 3, we have, by (30) and the symme- 
trical equations for G and H, 



+ dx 
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Gr=(m-n)h. & +nUr~^-iy+-^-\, 

Hr=(m-n)>.*+n{(r|-l) y +-^}. 



J 



♦ *.*,, (34) 



26. It is to be remarked that these equations express such functions of (#, y, z\ the 
coordinates of any point P of the solid, that F.^ G-*>, H.a> are the three components of 
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the force transmitted across an infinitely small area a perpendicular to OP, while, for any 
point of either the outer or the inner bounding spherical surface, Fa/, G#, Ha> are the 
three components of the force applied to an infinitely small element a of this surface. 

27. To reduce the surface-equations of condition derived from these expressions to 
harmonic equations, let us consider homogeneous terms of degree i of the complete 
solution, which we shall denote by a 4 , j3 i5 y i9 and let Vi*> Ki+i denote the correspond- 
ing terms of the other functions. Thus we have 



Gr=X< (m— n%__ ] y-\-n(i'—l)^ i ^ r n 



dx 
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H.r=%< 



dy 



(m^n)\^ l z+n(i—l)y i +n'-^ L 



(35) 
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28. The second of the three terms of order i in these equations, when the general 
solution of § 13 is used, become at the boundary each explicitly the sum of two surface 
harmonics of orders i and i— 2 respectively. To bring the other parts of the expressions 
to similar forms, it is convenient that we should first express £ i+1 in terms of the general 
solution (12) of § 13, by selecting the terms of algebraic degree i. Thus we have 
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and symmetrical expressions for ft and y i? from which we find 

Hence, by the proper formulae [see (42) below] for reduction to harmonics. 
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and (as before assumed in § 12) 

, dtii , dvi . dwi 
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Also, by (11) of § 12, or directly from (36) by differentiation* we have 
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Substituting these expressions for $ t _ l9 a { , and ^ + , in (35), we find 



* When t—1 is positive, <^_ai will express the same function as Y*_i of § 9 above. The suffixes now intro- 
duced have reference solely to the algebraic degree, positive or negative, of the functions, Whether harmonic or 
not, to the symbols for Which they are applied. 
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._*/. 



n(2i— 1) [(m— 2n)i 4- 2m + n\ 



n[2i(i— l)m— (2i— «l)n] 2 ^-i 



r 



Fr—2^— 1>,+ (2i + i)[( TO + 2»)i-»i-n] ^- 1 " fc (2i + l)[(m + 27i)i-m-7i]' ^ 

This is reduced to the required harmonic form by the obviously proper formula 



2i+l &r J* W 



^-i — 2 i-l* 



^ 



#,-_, - l+1 d(^ I r-«*+ , )l 



dx 



—r 



dx 



(42) 



Thus, and dealing similarly with the expressions for Qr and H>, we have, finally, 



Fr =nZUi-l)u { -2(»-2)M i r*%i-E,r' 



, +1 rf(fr-.T-" +I ) 



„2t+ 



ef<p 



i+i 



dx 

'■ ' ' ' "" ■ "«l 



d# 



2z +1 afc 



•p.jM-i 'ifrfo-i*-*") !_ fell 

' (te 2t'+l dy y 



Gr=»s{(»-1>, -2(*-2)M/ 
Hr=^(<-l)«>-2(*-2)M^^ %>} 



(43) 



where [as above, (16) of § 16.] 



1 



and now further 



* 2 {m J r2n)i— m— n 

(m--2n)i+2m + n 



JLJIi* ■■■■■• 



(2i H-l) [(m + 2n)i— m— n\ 



Y. 



. • • • • • 



(44) 



29. To express the surface conditions by harmonic equations, let us suppose the 
superficial values of F, G, H to be given as follows : 

►when r=a, 



and 



G=SB 



/ >, 



f=2a-, 

G=^B-, 

— ~4^i ? 



when r=&', 



(45) 



where A B i5 C f , AJ, B^, CJ denote surface harmonics of order i Now the terms of 
algebraic degree % exhibited in the preceding expressions (43) for Fr, Gr, Hr, become, 
at either of the concentric spherical surfaces, sums of surface harmonics of orders % and 
f— 2, when % is positive, and of orders — i— 1 and — i— 3 when £ is negative. Hence, 
selecting all the terms which lead to surface harmonics of order % and equating to the 
proper terms of the data (45), we have 



n 
r 



1 dx - Ej -i~i / ^ 2i+l 



dx 



<fy-A 



A* when r~#, "1 

*■=< * '• (46) 

A[ when r=a\\ 



and symmetrical equations relative to y and z. 

30 These equations might be dealt with exactly as formerly with the equations (15) 
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of § 15. But the following order of proceeding is more convenient. Commencing with 

the first of the surface equations (46), multiplying it by I- ) , attending to the degree of 

each term, and taking advantage of the principle that, if $ he any homogeneous function 
of #, y> z, of degree t, the function of angular coordinates, or of the ratios w:y:z, 

^ for any value of r, 



which it becomes at the spherical surface r=a, is the same as 
we have 




n 



(i_l) Mi ._(,- +2 ) I ^^.-M^^-f +2(i+l)M_ i+l a>('- 



2i+l 



dx 

d<p. 



■^ dx — -^-i-i« dx ~ 2i+lL<fe~Vv & -I 



=^(0'' ( 47 ) 



where the second member, and each term of the first member, is now a homogeneous 
function of degree % of #, y, # (being in fact a solid spherical harmonic of degree and 

order i). Taking -r- of this, and j- and -t- of the two symmetrical equations, adding, 

taking into account equations (38) and (39), and taking advantage of the equation 
V 2 V=0 for the solid harmonic functions concerned, we have 

I (48) 

1 fd(A^) , d{B^) . didriy ' v ' 



a { \ dx 



+ 



dy 



mmLmm 



dz 



)■ 



d 



Again, multiplying (47) by arV~ 2i ~V and taking ^ 2i+3 ^ of the result, dealing similarly 
with the two symmetrical equations, and adding, we have 



n 
— < 
a 



f / r \2i+3 ^ 

2m-> i+1 -[*+2-(2i+l)(i+l)E_ i _ 1 ]Q 4_ s _ a +2^+l)(2i+l)M i+ ^ i+1 J 



,2i+3 



f^Air-*- 1 ) , d(B*r-<->) . rf(<V-*-») 






+ 



} 



(49) 



dx ~ r dy ^ dz 
Changing i into i-— 2 in this equation, we have 

^2(i-2K^. 1 -P-(2i-3Xi-l)E_ i+ jQ^V-i+2(«-2)(i-lX2i-3)M^_ 1 } 



r^- 1 f djAj^r-^ 1 ) djBj^r-^ 1 ) d(Cj- 2 r- t+i ) 
a' \ dx ' dy ' dz 



I 



(50) 



Precisely similar equations, derived from the inner surface condition of the shell, are 
obtained by changing a, A, B, C into d, A', B', C\ We thus have (48), (50), and the 
two corresponding equations for the inner surface, in all four equations, to determine 
the four unknown functions v^, 4>-i> ft-i? <?-*> * n terms of the data which appear in 
the second members. The equations being simple algebraic equations, we may regard 
these four functions as explicitly determined. In other words, we may suppose <ft and 
$i known for every positive or negative integral value of L Then equation (47), the 
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two equations symmetrical with it, and the others got by changing A, a, &a into A', d<> 
&c., give Un v t , w ( explicitly in terms of known functions, and the expressions (36) 
for a i9 j3 t -, yi complete the solution of the problem. 

31. The solution for the infinite plane plate is of course included in the general solu- 
tion for the spherical shell, as remarked above for the case in which surface displace- 
ments, not surface forces, were given ; but, as in that case, it will be simpler and prac- 
tically easier to work out the problem ah initio, taking advantage of the appropriate 
Fourier forms. The relative ease of the independent investigation is indeed still greater 
in the case in which the surface forces are given than in the other case, since the general 
expressions for the surface forces assume simple forms when the surface is plane, and 
require no such transformation as that which we have found necessary, and which has 
constituted the special difficulty of the problem, when the surface was spherical. The 
problem of the plane plate presents many questions of remarkable interest and practical 
importance; and although the object and limits of the present paper preclude any 
detailed investigation of special cases, we may make a short digression to work out the 
general solution. 

32. Let the origin of coordinates be taken in one side of the plate and the axis OX 
perpendicular to it. Then, according to the general expressions (25) of § 21, the three 
components of the force per unit of area, in or parallel to either side of the plate, are 
respectively 

dy 

dz 



parallel to OX, P=(m+%) -r-+(m— ^)(;r+" 
parallel to OY, U=»(|+g 
paraUeltoOZ.Tm+g) 



(51) 



The surface condition to be fulfilled is that each of these functions shall have an arbi- 
trarily given value at every point of each infinite plane side of the plate. 

33. From the indications of § 6 above, it is easily seen that the following assump- 
tions are correct for a general solution of the equations of internal equilibrium, and con- 
venient for the application at present proposed, 

d<p 
1 ax 

d<p 

where % v, w, and <p denote functions of (#, y, z) which each fulfil the equation V 2 V=0. 
From these, by differentiation, and by taking y 2 <p=0 into account, we have 

dot x dp | dy du dv dw df 

dw*dy~*~dz dw*dy* dz * doe 5 
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or 



p— >ji.i. — 



if 



dx 



. ^du do dw 
y~~~da;*dy*dz ' 



(52) 



and I be used with the same signification as above (§ 2). Also, by differentiation and 
application of the equation, V«=0, V |=0, we find 

dx 1 ' r dxdy ' ' dxdz 

Hence, to satisfy the general equations of internal equilibrium (3) of § 3, we must 
have 

dtp m 



dx' 



m+2n 



*. 



Hence the general solution becomes 



mx 



ccz=u — 



$=v — 



m-\-2n 



mx 



dvfydi 



<y=w 



m-\-2n dy 
mx dfydx 



m + 2n dz 



y • 



• •••• ••• i^jtji 



j 



where u, % w are any functions whatever which satisfy the general equation V 2 V=0, 
and \f/ is given by (52) ; and where, further, it must be understood that fydw must be 
so assigned as to satisfy the equation V 2 V=0, which \|/ itself satisfies by virtue of (52). 
34. The general form of the solution of V 2 V=0, convenient for the present applica- 
tion, is clearly 

where jp, s, t are three constants subject to the equation 

2 2 I -42 

JO -— O ""J"" V • 

If now we suppose, as a particular case, the surface condition to be that 



and 



P= A sin (sy) sin (fe), , 
U=B cos (sy) sin (tz\ >when x 
T=C sin (sy) cos (fo), 



=0. 



HDCCCLXIII. 



P=A'sin (sy) sin (fe), 

U=B'cos (sy) sin (fe), 

T=C ; sin (sy) cos (fe), 

4 M 



when %=a, 



. . . (54) 
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where A, B, C, A', B', G are six given constants, we must clearly have 

u=(fi-**+f 1 ?*) ffln'(«yj sin (fe), 

v=(gs- p *-\-g , e? x ) cos (sy) sin (fe), 

w=(kr px -\-hfe?*) sin (sy) cos (fe), 



(55) 



where /, g, % f, g\ 72 are six constants to be determined by six linear equations 
obtained directly from (54), (51), (53), (52), (55). Btit, by proper interchanges of 
sines and cosines, we have in (54) a representation of the general terms of the series 
or of the definite integrals, representing, according to Fourier's principles, the six 
arbitrary functions, whether periodic or non-periodic, by which P ? U, T are given over 
each of the two infinite plane sides. Hence the solution thus indicated is complete. 

35. To complete the theory of the equilibrium of an elastic spheroidal shell, we 
must now suppose every point of the solid substance to be urged by a given force. 
The problem thus presented will be reduced to that already solved, by the following 
simple investigation. 

36. Let X, Y, Z be the components of the force per unit of volume on the substance 
at any point #, y, z. (That is to say, let §X, #Y, qZ be the three components of the 
actual force on a volume q, infinitely small in all its dimensions, enclosing the point 
'(#, #, z). Not to unnecessarily limit the problem, we must suppose X, Y, Z to be each 
an absolutely arbitrary function of #, y, z. 

37. When we remember that #, y, z are the coordinates of the undisturbed position 
of any point of the substance, and differ by the infinitely small quantities a, |8, y 
from the actual coordinates of the same point of the substance in the body disturbed 
by the applied forces, we perceive that JLd%-\-Ydy-$-Zdz need not be the differential 
of a function of three independent variables. It actually will not be a complete 
differential if the case be that of the interior kinetic equilibrium of a rigid body 
starting from rest under the influence of given constant forces applied to its surface, 
and having for their resultant a couple in a plane perpendicular to a principal axis. 
Nor will ~X.dw-\-Ydy-\-Zdz be a complete differential in the interior of a steel bar- 
magnet held at rest under the influence of an electric current directed through one 
half of its length, as we perceive when we consider Faraday's beautiful experiment 
showing rotation to supervene in this case when the magnet is freed from all mecha- 
nical constraint. 

38. The equations of elastic equilibrium are of course now 



riV 2 a 4-m -r— — X, 

ax 
dy 



» •••••«•«• I <DO J 
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Let 7S, g, g denote some three particular solutions of the equations 



W=-X, 



V 2 * 



Y, 



» • • • • • • • • \ J 



vv=-z. 

These, sr, g 9 tr, we may regard as known functions* being derivable from X, Y, Z by known 
methods (Thomson and Tait's ' Natural Philosophy,' chap. vi.). Then, if we assume 



1ST 

n 



/> 



£_ 






r--=y, 



(58) 



/ 



and 



dx * dy*dz ' 5 ' 
the equations (56) of interior equilibrium become 

' ' ax n dx 

T-ron . d$. m d£ 

"X^+ m Ty = -nTy> 

— -„ , d§. m d£ 

where £ is a known function given by the equation 

* dx*dy*~dz 



(59) 



• • • • • « • • 



(61) 



Now, as we verify in a moment by differentiation, equations (59) and (60) are satis- 
fied by - m j® V 



a 



1 n{m~\-n) dx 
•—m d$ 



ft= 



n(m + n) dp 

—m d§ 

*i~n{m*\~n)dz , 



(62) 



if & is some particular solution of 

V 2 ^=| (63) 

Hence (58), (57), (62), (63), (61) express a particular solution of (56). 

39. We conclude that the general solution of (56) may be expressed thus: — 



1/ 

n\ m-\-n ax 



m d§\ . ^ 



P=l(i 



Y=Z\ ff ~ 



m d$ 



n 



m-t-n dy 

m dh 
m-\-n dz 

4M 2 



) +v> 

' ) 



(64) 
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where 



w =v- 2 x, 

<r =V" 2 Z , 



/dm dg d 



§=V~ 2 {~+ 



\dx 






dy'dz 



(65) 



according to an abbreviated notation, which explains itself sufficiently ; and 'a, ; /3 5 'y 
denote a general solution of the equations 



nV u~\-m-r- 



d ld ] a -d'P . d 



+-^54 



dx \dx* dy 
i-fof/i . d /d'ct t d f 8 



+ 



dz 
d 



v-,or i d /d ! ct 



h^4 



^2/ 



£)=* 



(66) 



/ 



40. This solution is applicable of course to an elastic body of any shape. It enables 
us to determine the displacement of every point of it when any given force is applied 
to every point of its interior, and either displacements or forces are given over the 
whole surface, if we can solve the general problem for the same shape of body with 
arbitrary superficial data, but no force on the interior parts. For 'a, 'j3, 'y are deter- 
mined by the solution of this problem, to be worked out with the given arbitrary super- 
ficial functions modified by the subtraction from them of terms due to the parts of a, j3, y 
which are explicitly shown in terms of data by equations (64) and (65). 

41. Hence the problem of § 35 is completely solved, — whether we have displacements 
given over each of the two concentric spherical bounding surfaces, when the solution of 
§§ 14-18 determines 'a, '/3, 'y ; or forces given over the boundary, when the solution of 
§§ 26-30 is available. In the former case the superficial values of the functions 



1 / m d$\ 

ny m+ndyj' 



m + ndx 
d$ 



1/ m d$\ 

n\ m+ndzV 



known from equations (65), must be subtracted from the arbitrary functions given as 
the superficial values of a, /3, y, and the residues, expressed in surface-harmonic series 
by the known method, will be the harmonic expressions for the superficial values 
of 'a, '/3, 'y. In the latter case, we must first substitute those known functions 

_/ OT \ & c#) instead of a, /3, y respectively in (34), and the values of Fr, Gr, H> 

thus found must be subtracted from the given arbitrary functions representing the true 
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superficial values of Fr, Gr, Hn The remainders, which we may denote by f Fr, 'Gr, 'H>, 
must then be reduced to harmonic series, as in (45), and used according to the investi- 
gation of § 30, to determine 'a, '/3, 'y. 

42. The general solution (64) and the expression just indicated for the terms to be 
subtracted from the data so as to find Tr, ! Gr , 'H> , becomes much simplified when, as 
in some of the most important practical applications, ILdtt-^Ydy-j-Zdz is a complete 
differential. Thus let 



^_dW V _^W nr_dW 

da? ay dz 



(67) 



W denoting any function of x, y, z. Then, assuming, as we may do according to (65), 

d 



-=&v- 2W ' 



d 



e=T y 



v- 2 w, 



d 



we have by differentiating, &c. ; 



and therefore 



Hence the solution (64) becomes 



d^sr . dg . dcr _^_ 
doe "•" dy* dz ' 

$=v~ 2 w. 



(68) 



1 d$ , i 

a = — . — —. ^'^ 



(*■ 



m-\-n dx 
1 d$ 



m + ndy 
1 d& 



+% 



* m-\-n dz* ' * 



• • 



. . (69) 



From this we find 



and (§ 25) 






* 



m-\-n 



w+% 






> 



and 



Hence, by (34), 



But 



/> d f oc | d f 3 , rf'y' 
dx dy dz 



• • (70) 



Fr. 



m+n 



(m-n)Wx+n[(rz-l)£+&r %]$}+& 



r dr~ \dr~T~ L )dx 



d_ d 
dx 
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Thus for Fr, and the symmetrical expressions, we have 

1 f - Ad *X 



Fr 



Gr 






Hr= \ (m—ri)Wz +2nr 



m + n 



dr dz 



+'Hr. 



) 



. . . (71) 



43. These expressions become further simplified if W is a homogeneous function of 
any positive or negative integral or fractional order £+1? in which case we shall denote 

it by W i+1 . For 9- will be a homogeneous function of order i+3, and ^ of order i+2. 
Hence 



Hence the preceding become 



4s*=(*+2) 



dx 



Fr=^(m-n)W i+1 a?+2n(*+2)§}+'Fr, > 



Gr 



Kr 



m-\~n 

1 

m-\-n 

1 

1 ' ""■■■ « 

m-\-n 



(m— n)W i+l y+2n(i+2) 



dy 
£»1 



+ f Gr, 



(w— n)W i+1 z-\-2n{i-\-2)-^ 



+'Hr. 



.....( I u ) 



) 



44. These expressions are the more readily reduced to the harmonic forms proper for 
working out the solution, if the interior force potential, W m , is itself a harmonic 
function. We then have (§ 10) 

1 



&: 



and 



which give 



2(2^ + 5) 



*.2 W ** — /Vw -L.l r 2dWi+1> \ 

r VV * +1 > dx~ 2i + 5\ Wyy ^ 1 ^ r 2 r dx ) 



W i+1 x-. 



' dW i+1 



2i + £\ T dx 



nn 



2i+5 



d(W i+1 r-^ 
dx 



(73) 



Vr 1 f «i+(»+l)n rfW l+1 w(2i+5) -n i+5 d^^^^ ) 






(74) 



and symmetrical expressions for Gr and Hn Here the terms to be subtracted from the 
arbitrary functions given to represent the superficial values of Fr, Gr, and Hr are each 
explicitly expressed in sums of two surface harmonics of orders i or — i— 1, and i-\-2 or 
—i— 3 respectively, viz., in each case, that one of the two numbers which is not negative. 
45. When the shell is in equilibrium under the influence of the forces acting on it 
through its interior, without any application of force to its surface, we must have 

Fr=(V 

Gr = 0, hvhen r=a and when r=a f (75) 

Hr=0, 
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Hence, for the case in which W is a spherical harmonic, the preceding equations 
give the proper harmonic expressions for 'Fr, ; Gr, 'B> at the outer and inner bound- 
ing surfaces, for determining ; a, f j3, 'y by the method of §§ 28-30. Thus, using all 
the same notations, with the exception of 'a, '|3, 'y» 'F, '^ 'H, instead of a, j3, y, 
F, G, H, and, for the present, supposing i+1 to be positive*, we have the complete 
harmonic expressions of 'F, f G, f H, each in two terms, of orders i and e+2 respect- 
ively. Hence the A, A', &c, of (45) are given by the following equations ; — 



A* _ A! ____ »t+(i + l)fl g-i jWf+i 

da? 



."+ 1 (2i + 3)(m+n) 



B* B, 



t + 1 



wi +(* + !)» y -, 5 



a 



^4 , 



a 



' i+1 (2i + 3)(w + ») 



% 



C« ______ w+ (»+!)» -<dW< +1 



a 1 



/<+i -(2i + 3)(j» + n) 



*fe? 



A<+a_ AJ +8 _ (2i+5)ro— w yi+s ^Wi+.r-*- 3 ) , 
^+i a «+i (2i+3)(2i+5)(»i+») & 

B i+s _B; +a _ (2* + 5)w-n .^+8^(W <+ ir- af - s ) 



^•+i a /«+i (2t + 3)(2i + 5)(m+n) 



/M» 



«fy 



c- 



f+2 



C- 



i+2, 



(2i + 5)i»-» ^ +8 <?(W <+1 r-*-*) 



fl<+ 1 a?*+* (2i+3)(2i + 5)(m + w) 



r 



afe 



/ 



(76) 



46. The functions derived from A,., B i5 Q, &c, which are required for formulae (48) 
and (49), are therefore as follows : — 



dx 



+ 



dy 



.d(Crf 
"+" dz 



diAtr-i- 1 ) , d(B___) i((V-'-') 



U/ffi 



dy 



dz 



d(A <+ _r < + 2 ) J(B< +a H +a ) ^((W i+2 ) 

+ ^, "T 



fifo? 



% 



% 



o, 



(t + l)(2t + l)[w+(t + l)w] fl_^ w 

(2i + 3)(m+») r*+* W ' +1 ' 

(f + 2)[(2i + 5)ro—n] i+1 w 



V. 



<*(A, +£ r-«-») J(B i+2 r-^ 3 ) + ^(C, +2 r-^) 



dfo 



dy 



dz 



=o, 



/frfrv 



/ 



with the corresponding expressions relative to AJ, B^, C£, &c, obtained simply by 
changing a into a!. 

Hence by (48) and (50), and the two corresponding equations for the inner surface, 
we infer that each of the four functions \}/ i _ 1 , ^, <p t „ 15 <p_ t . vanishes. By the same 
equations, with i changed into $+2, we obtain expressions, all of one harmonic form, 
direct or reciprocal, as follows, for the four functions of order i+1 :— 



* As we shall not in the present paper consider particularly any case of a shell influenced by centres of 
force in the hollow space within it, which alone could give a potential W* + i of negative degree, we need not 
write any of the expressions in forms convenient for making i+1 negative. 
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(78) 



K i+1 , KJ +1 , L i+1 , L) +1 , which need not be here explicitly expressed, being four constants 
obtained from the solution of four simple algebraic equations. Lastly, by the four 
equations with (£+4) instead of i 9 we find that \J/ i+3 , %J/_j_ 4 , <p i+3 , <p_i_ 4 all vanish. Using 
these results for *p and <p in (47), we see that each of the functions u must be a harmonic 

dW- d(W- r" 2 *'"" 3 ) 
congruent with either / +1 or — — ^ -• Hence, by using (78) in (38) and (39) we 



find 



u t 



__. —L i+1 dW i+l 
(2i+l)(i+l) <fo 



flfo 



-K 



*+i 



*+2 



(2t + 5)(» + 2) 



/y»' 



»+, rf(W <+1 r-*-») 
dx 



- f " 1 (2t+l)(t+l) & ' I 

^L; +1 <Z(W i+1 r- 2 - 3 ) . 



w„i_ 9 = 



(2t + 5)(t + 2) 



cfc 



(79) 



and symmetrical expressions for v and w. Finally, using these expressions, (79) and 
(78), in (36), and the result in (69) with (73), we arrive at an explicit solution of the 
problem in the following remarkably simple form : — 



„_« «*W, +1 , - rf(W <+ ,!-«-») 



rfy 



rfy 



^-e «*W« + . . tf rf(W <+1 r-»-») 



• . 



. . (80) 



where 



C 



i+1 






(t+l)(2i+l) ^1 



2(2* + 3)(m + n) 



— 2 m 



K 



i+i 



(m + 2ft)z + m + 3tt 



W, 



__ K,+,i* + » + L 



»2i+5 



(f + 2)(2i + 5) 



+ 



«if 



K'. 



i+1 



+ 3)(2i + 5)(m + ft) ' 2 (m + 2w)i+2m + 3ft 



}• 



(81) 



47. In conclusion, let us consider the case of a solid sphere. For this we have 

%|/_<_ a =0, and <p„ f _ 2 =0, 

as we see at once from the character of the problem, or as we find by putting a f =0 in 
the four equations by which in § 46 we have seen that K l+1 , KJ +n L i+1 , L* +1 are to be 
determined. Then, by (48), with i changed into £+2, and by (49), we find 



f _ (i+2)[(ro + 2ft)i + m + 3w][ira(2i+5)-- n] w 
Yi+i— "~(2i + 3)w{/w[2(i + 2) 2 +l]-w(2i + 3)}(m + ^) i+1 ' 

^ o (»>l)»(2»>l)[m(i>3)-n] 

9i + 1 "-^2^{m[2(i + 2) 2 +l]~w(2i + 3)} . vv «+»- 



• (82) 
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The coefficients of W <+1 in these expressions are the values which we must take for 
Kf +1 and L,. +1 respectively in (81); and therefore, after reductions which show (m+^) 
as a factor of the numerator of each fraction in which it appears at first as a factor of 
the denominator, we have 



©,•4-: = 



(i + l)[>(i + 3)-n]a< 



'i+ 






[(i + 2)(2i + 5)wi-(2i + 3)w]r s 



®m=; 



2n{m[2(i + 2)* + l]~n(2i + 3)}^27i{2i + 3){m[2{i + 2)*+ll--n{2i + 3)} 
(z+l)mr 2 * +5 



. (83) 



*w, +I 



n{2i + 3) {m[2{i+ 2) 2 + 1] — n(2i + 3) } 

These, substituted in (80), give expressions for a, /3, y which constitute a complete and 
explicit solution of the problem. 

It is easy to verify this result, by testing that (56) (with -~X=^p^ , &c.) is satis- 
fied for every point of the solid, and that equations (34) give F=0, G=0, H=0 at the 
bounding surface, r=a. 

48. The case of i=l is, as we shall immediately see, of high importance. For it the 
preceding expressions, (83) and (80), become 

__ — 10(4m— w)a 2 -f (21m— bn)i a ] 



«,= 



10n(l9m-~5ri) 



Ci = 



4mr 7 



(%•< 



10n(l9m— 5n) 

(XJLf G/Qff 



V 



(84) 



M 



P~^ 2 dy ^ 2 dy ' 
' dz as 

49. As an example of the application of §§45-48, let us suppose a spherical shell or 
solid sphere to be equilibrated under the influence of masses collected in two fixed 
external points*, and each attracting according to the inverse square of its distance. 
Let the two masses M, M' be in the 

axis OX; and, P being the point whose M f 

*■—' i i n . m . 1 i ■ i ii i i 

coordinates are #, y, 2, let PM=D, 
PM'=D f . Let also OM=c, OM!=cf. 
Then, if m, m' denote the two masses, for equilibrium we must have 

m m! 

c 2 t J< * 

* If our limits permitted, a highly interesting example might be made of the case of a shell under the 
influence of a single attracting point in the hollow space within it. The effect will clearly be to keep the 
whole shell sensibly in equilibrium eyen if the attracting point is exeentrie; and under stress even if the 
attracting point is in the centre. 
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m . m' 



The potential at P, due to the two masses, will be ]j+w? or, according to the notation 
of § 42, with, besides, w taken to denote the mass of unit volume of the elastic solid, 

The known forms in the elementary theory of spherical harmonics give immediately 
the development of this in a converging infinite series of solid harmonic terms. We 
have only then to apply the solution of §§ 45, 46 to each term, to obtain a series 
expressing the required solution. 

50. We may work out this result explicitly for the case in which both masses are 
very distant; and for simplicity we shall suppose one of them infinitely more distant 
than the other ; that is to say, we shall suppose it to exercise merely a constant balancing 
force on the substance of the shell. We shall then have precisely the same bodily dis- 
turbing force as that which the earth experiences from the moon alone, or from the sun 
alone. 

51. Referring to the diagram and notation of § 49, we have 



D — c 
x y 



1+-+ ^ ; 



if we neglect higher powers of -, -, - than the square ; and 

Lt Lt (s 

neglecting all higher powers of £ , y - , -. Hence, taking account of the relation ~=^ 
required for equilibrium, we have, for the disturbance potential, 



c? 2 c' 2 



ui 



an irrelevant constant being omitted from the expression which § 49 would give. This 
being a* harmonic of the second degree, we may use it for W i+ „ putting i=l in the 
formulae of § 47, and thus solve the problem of finding the deformation of a homo- 
geneous spherical shell under the influence of a distant attracting mass and a uniform 
balancing force. I hope, in a future communication to the Eoyal Society, to show the 
application of this result to the case of the lunar and solar influence on a body such as 
the earth is assumed to be by many geologists — that is to say, a solid crust, constituting 
a spheroidal shell, of some thickness less than 100 miles, with its interior filled with 
liquid. The untenability of this hypothesis is, however, sufficiently demonstrated by the 
considerations adduced in a previous communication (« On the Eigidity of the Earth," 
read May 8, 1862), in which the following explicit solution of the problem for a homo- 
geneous solid sphere only is used. 
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52. Using the expression of § 51 for W 2 , we have 









m 



+ 5T3W, 






m 

+ -3ZW1 



► , (85) 



^(W 2 r-^) __ m (2^~^~-^% g£J __ q m (2g*- &«-£*«)* 

These formulae being substituted for the differential coefficients which appear in (84), 

we have algebraic expressions for the displacement of any point of the solid. 

The condition of the body being symmetrical about the axis of #, we may conveniently 

assume . ^ 

y=ycos£>, z=ysm<p, 

so that we shall have (as we see by the preceding expressions) 

/3=^cos<p, 
y=p sin <p, 

and (Jb will denote the component displacement perpendicular to OX. If, further, we 

assume M M ^ aA 

#=r cos#, 

y=rsin#, 
the expressions (84) for the component displacements, with (85) used in them, give 



c 3 



„ m 



3 % 

2 r A 



-2re 2 +~^ 2 (5cos 2 0--3) 



cos 0, 



r ® 2 +|§(5cos 2 ^-.l)|sin^. 



2 r ■ 



(86) 



The values given in (84) for ® 2 and $ 2 are to be used for any internal point, at a distance 
r from the centre, in these equations (86), and thus we have the simplest possible ex- 
pression for the required displacement of any point of the solid. 

53. If we resolve the displacement along and perpendicular to the radius, and con- 
sider only the radial component, we see that the series of concentric spherical surfaces 
of the undisturbed globe become spheroids of revolution in the distorted body. The 
elongation of the axial radius, obtained by putting 0=0 and taking the value of a, is 
double the shortening of the equatorial radius, obtained by putting 0=^t and taking 
the value of {h\ which we might have inferred from the fact, shown by the general 
equations (80) above, that there can be no alteration of volume on the whole within any 
one of these surfaces. The expression for the excess of the axial above the equatorial 



radius is 



4n2 
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which, if we substitute for <& 2 and (&' 2 their values by (84), becomes 

~ m 2 (4m— w)« 2 —(3m— n)r z 
W l? I0n{l9m--5n) T ' 

If in this we take r=a, and m=co , it becomes y^Trf 3 , which is the result used in 

§ 34 of the paper on the Eigidity of the Earth, preceding the present in the Transactions. 
54. In the case of a f =0, the result of § 18 takes the extremely simple form 



a — XiAti a J + 2fl «[ n (2,-_i)+m(i--l)] 



MB + 



m(a?-~ r 2 ) 



r 



7 



: M Q ( ~ 



where 



2# £ [ft(2i— l)-fwi(i— 1)] dy 

m(^-r 2 ) d@ I _ 1 

+ 2a i [n(2i— l)+m(« — 1)] & 



h 



(87) 



£) 
^*-i 



^A^) 



W. + <*(<V) 



dx l dy % dz 

This expresses the displacement at any point within a solid sphere of radius #, when its 
surface is displaced in a given manner (SA|, 2B t , SQ). And merely by making i 
negative we have, in the same formula, the solution of the same problem for an infinite 
solid with a hollow spherical space every point of the surface of which is displaced to a 
given distance in a given direction. These solutions are obtained directly, with great 
ease, by the method of §§ 6-15, or are easily proved by direct verification, without any 
of the intricacy of analysis inevitable when, as in the general investigations with which 
we commenced, a shell bounded by two concentric spherical surfaces is the subject. 



[Added since the reading of the Paper.] 

§§ 55 to 58. Oscillations of a Liquid Sphere. 

55. Let V be the gravitation potential at any point P(#, y, z), and h the height of the 
surface (or radial component of its displacement) from the mean spherical surface at a 
point E in the radius through P. Then, if 

be the expression for h in terms of spherical surface harmonic functions of the position 
of E, and if \h be the attraction on the unit of mass exercised by a particle equal in mass 
to the unit bulk of the liquid, we have, by the known methods for finding the attractions 
of bodies infinitely nearly spherical (Thomson and Tait's 'Natural Philosophy,' chap. vL), 



V=4?rap \a— -^—+2 



r 
a 



and 



V Si 



when r<a< 



33 



33 



r>a. 



r=a. 



(89) 
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In these 

4<z>pa=ty, . (90) 

if g denote the force of gravity at the surface, due to the mean sphere, of radius a. 

56. Now for infinitely small motions the ordinary kinetic equations give 

~dx~~$\dt dx) ; dy~~£\dt dy ) > dz~2\dt~~ dz )' ' ' ( 91 >> 

where g is the mass per unit of volume ; u,v 9 w the component velocities through the fixed 
point P at time t; and p the fluid pressure. Hence, possible non-periodic motions being 
omitted, udfi+vdy+wdz is a complete differential ; and, denoting it by <p, d we have 

C-*=»(4-v) '. (92) 

57. To find the surface conditions, — first, since the pressure has a constant value, n, 

at the free surface, 

_p=^A+n when r=a 9 (93) 

the variations of gravity depending on the variations of figure being of course neglected 

in the infinitely small term ggh. And, since jt is the radial component of the velocity 

at E, we have, when r=a, 

£L it l.1. fyji.L it-J^L . (94) 

r dx~^~ r dy~*" r dz dt 



Now since, the fluid being incompressible, V <p=0, <p may be expanded in a series of 
solid harmonic functions ; let 

*=**<($> • (96) 

where 2 , <E> 2 , . . . are surface harmonics. Hence, as the successive terms are homoge- 
neous functions of the coordinates (#, y, z), of degrees 1, 2, &c., 

r dx*r dy*r dz r \ a / * * * * * v 17 ^/ 

and therefore, by (88) and (94), 

*- = ;*« • ( 97 > 

58. Eliminating jp between (92) with r=a and (93), substituting for V by (89) and 

(90), differentiating, substituting for ~-^ by (97), and comparing harmonic terms of 
order i, we have 



of which the integral is 



^ ( =Acos 



{*v/f«(i-2iTl)- E } • W 



Here A is a surface spherical harmonic function of the coordinates of E expressing the 
maximum value of <E> <3 and E is the epoch (Thomson and Tait, § 53) of the simple 
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harmonic function of the time which we find to represent <£>,. Using this solution in 
(97) and (88), we see that if the surface be normally displaced according to a spherical 
harmonic of order i, and left to itself, the resulting motion gives rise to a simple 
harmonic variation of the normal displacement, having for period 



2i7f 



V 



a 2«+ 1 



, „, , (9i4-l)a -r 

that is, the period of a common pendulum of length vw 7 ~- • It is worthy of remark 

that the period of vibration thus calculated is the same for the same density of liquid, 
whatever be the dimensions of the globe. 

For the case of i=2, or an ellipsoidal deformation, the length of the isochronous 
pendulum becomes f #, or one and a quarter times the earth's radius, for a homogeneous 
liquid globe of the same mass and diameter as the earth ; and therefore for this case, or 
for any homogeneous liquid globe of about 5 J times the density of water, the half-period 
is 47 m 12 s , which is the result stated in the paper " On the Eigidity of the Earth" (§ 3), 
preceding the present in the Transactions. 

Appefdix, §§ 59-71. — General Theory of the Equilibrium of an Elastic Solid. 

59. Let a solid composed of matter fulfilling no condition of isotropy in any part, and 
not homogeneous from part to part, be given of any shape, unstrained, and let every 
point of its surface be altered in position to a given distance in a given direction. It is 
required to find the displacement of every point of its substance in equilibrium. Let 
#, y, z be the coordinates of any particle, P, of the substance in its undisturbed position, 
and w+a 3 y +/3, z+y its coordinates when displaced in the manner specified ; that is to say, 
let &, j3, y be the components of the required displacement. Then, if for brevity we put 



A -(£+ 1 ) +(s) + \£) 






C 



a-. 



*y. 






'dt 



OL 



dz 



dy 



dx) 

dy x 2 



+ (^+i 



\dz 



rh, Jv ~T \dy ~T jdz~^~ dy \dz * 



I— 



den, da 
dy dz 

da / den, 
d \doc 



«=(b+i) 



da dfi (dfi 



dz* ) dx J 



dy dy . 



dy l dw\dy l J ' doc dy 



(100) 



these six quantities A, B, C, a, b, <?, as is known * thoroughly determine the strain expe- 
rienced by the substance infinitely near the particle P (irrespectively of any rotation it 
may experience), in the following manner : — 

* Thomsok and Tait's c Natural Philosophy/ § 190 (e) and § 181 (5). 
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60. Let i, tj 9 £ be the undisturbed coordinates of a particle infinitely near P, rela- 
tively to axes through P parallel to those of #, y, z respectively ; and let § y , n i , £, be 
the coordinates, relative still to axes through P, when the solid is in its strained con- 
dition. Then 

^+^.+^=A?+B^+C^+2^+2^g+2cfo; .... (101) 

and therefore all particles which in the strained state lie on a spherical surface 

are 5 in the unstrained state, on the ellipsoidal surface, 

This, as is well known*, completely defines the homogeneous strain of the matter in 
the neighbourhood of P. 

61. Hence the thermo-dynamic principles by which, in a paper on the Thermo- 
elastic Properties of Matter in the first Number of the 4 Quarterly Mathematical Journal' 
(April 1855), Green's dynamical theory of elastic solids was demonstrated as part of the 
modern dynamical theory of heat, show that if w . dxdydz denote the work required to 
alter an infinitely small undisturbed volume, dxdydz, of the solid, into its disturbed con- 
dition, when its temperature is kept constant, we must have 

w=f(A 9 B y C 9 a,b 9 c), (102) 

where/ denotes a positive function of the six elements, which vanishes when A— -1, 
B— 1, C— 1, a, #, c each vanish. And if W denote the whole work required to produce 
the change actually experienced by the whole solid, we have 

W=SS$wdwdydz, . (103) 

where the triple integral is extended through the space occupied by the undisturbed 
solid. 

62. The position assumed by every particle in the interior of the solid will be such as 
to make this a minimum, subject to the condition that every particle of the surface takes 
the position given to it, this being the elementary condition of stable equilibrium. 
Hence, by the method of variation, 

lW=$$?>wdxdydz=0 (104) 

But, exhibiting only terms depending on Xa, we have 



to-k— (— + 1^4-— -+-- 
0W — \?dA.\dx + V + db dz^dc dy 

( O dwdot dw dot dw /dot - \ 



dlot 
dx 

dlot 
dy 



+ 



'dw dot dw dot dw /dot - \ 1 ddot 
dC dz^~ da dy*~ db\d%* ) J dz 



+&c. 
* Thomson and Tait's < Natural Philosophy/ §§ 155-165. 
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Hence, integrating by parts, and observing that &», ^3, ly vanish at the limiting surface, 
we have 

VN=-)))dxdydz (^-+^r+ir)5«+&c.k .... 



"1 \<fe ' dy' dz ) 



(105) 



dla diet dla 

J — -r- ? 
ft*/ 



where for brevity P, Q, E denote the factors of -^ ? ^5 ^ respectively, in the pre- 

ceding expression. In order that %W may vanish, the factors of &a, cS(3, cJy in the expres- 
sion now found for it must each vanish; and hence we have, as the equations of 
equilibrium, 



d f 9 dw /da. \ dw du dw dot 
d%[ dA\d%'* J*db dz*dcdy 

d (^dw da . diu da . dw /da 



JLJL-J 9 ZZ. ZZ.JL.ZH. HZ f , : 

dy\^dB dy* da dz* dc \d% 

d f rJtw dec , dw da . dw /da 



—j— JL 



(106) 






4 



J^A 






1. 



0, 



tfej^dC dz n ~ da dy~ r db \dx~ T ~' i ~J j 

iDCijm 43CC#« 

of which the second and third, not exhibited, may be written down merely by attending 
to the symmetry. 

63. From the property of w that it is necessarily positive when there is any strain, it 
follows that there must be some distribution of strain through the interior which shall 

make ^wdwdydz the least possible , subject to the prescribed surface condition, and there- 
fore that the solution of equations (106), subject to this condition, is possible. If, what- 
ever be the nature of the solid as to difference of elasticity in different directions, in any 
part, and as to heterogeneousness from part to part, and whatever be the extent of the 
change of form and dimensions to which it is subjected, there cannot be any internal 
configuration of unstable equilibrium, or consequently any but one of stable equilibrium, 
with the prescribed surface displacement and no disturbing force on the interior, then, 
besides being always positive, w must be such a function of A, B, &c. that there can be 
only one solution of the equations. This is obviously the case when the unstrained solid 
is homogeneous. 

64. It is easy to include, in a general investigation similar to the preceding, the 
effects of any force on the interior substance, such as we have considered particularly 
for a spherical shell, of homogeneous isotropic matter, in §§ 35-46 above. It is also 
easy to adapt the general investigation to superficial data of force^ instead of displace* 

65. Whatever be the general form of the function f for any part of the substance, 
since it is always positive it cannot change in sign when A— -1, B— 1, C— 1, a, #, c have 
their signs changed ; and therefore for infinitely small values of these quantities it must 
be a homogeneous quadratic function of them with constant coefficients. (And it may 
be useful to observe that for all values of the variables A, B, &c, it must therefore be 
expressible in the same form, with varying coefficients, each of which is always finite, for 
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(108) 



all values of the variables.) Thus, for infinitely small strains, we have Gbeen's theory 
of elastic solids, founded on a homogeneous quadratic function of the components of 
strain, expressing the work required to produce it. Putting 

A-l=2tf, B-l=2/; C~-l=2g, . (107) 

and denoting by ^(0, e\ ^{f,f\ . . . (#>/)» ...(£,«),... the coefficients, we have 

^= Ufa *y+(f> /)/*+& g)g 2 +(*, ay +(J, 5)5* +(c, c>? 2 } ^ 

+( e >f) e f+( e > g) e 9 +(*> «>« +(*> *>* +(^ e)«? 
+(/> ^+(/ ^+(/ 5 »l^+(/, ^ 
+(& a )g a +(g> %* +(g> c )g c 

+(#, V)db +(#, <?)#<? 
+(#> c)bc 

The twenty-one coefficients in this expression constitute the twenty-one coefficients of 
elasticity, which Green first showed to be proper and essential for a complete theory of 
the dynamics of an elastic solid subjected to infinitely small strains. 

66. When the strains are infinitely small, the products -^ ~~> -4f -A &c. are each 

J r dk dx db dz 

infinitely small, of the second order. We therefore omit them ; and then, attending to 
(107), we reduce (106) to 



d dw d dw d dw 



v 



dw de 

d dw 
dx dc 

d dw 
dx db 



d dw 



dz db 

d dw 
dz da 



d dw d dw 
' dy da ' dz dg 



0, 



0, 



0, 



(109) 



which are the equations of interior equilibrium. Attending to (108) we see that 
• • • are linear functions of e, f, g, a, b, c the components of strain. Writing out 



dw dw 

• • *. _____ 

de da 

one of them as an example, we have 

-£=(e 9 e)e+(e,f)f+fag)g+(e, a)a+(e, b)b+(e, c)c. 



. (110) 



And a, |3, y denoting, as before, the component displacements of any interior particle, P, 
from its undisturbed position (xy, z\ we have, by (107) and (100), 



dx J dy 



*=_■• 



„—A?> i dy t_.dy.dx 
dz dy dx dz 



d* dp 
dy~*dx 



(111) 



It is to be observed that the coefficients (#, e) (#,/), &c will be in general functions of 
(?, y, z\ but will be each constants when the unstrained solid is homogeneous. 
mdccclxiii. 4 o 
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67. It is now easy to prove directly, for the case of infinitely small strains, that the 
solution of the equations of interior equilibrium, whether for a heterogeneous or a homo- 
geneous solid, subject to the prescribed surface condition, is unique. For let a, j3, y be 
components of displacement fulfilling the equations, and let a', j3', y' denote any other 
functions of (#, y, z) having the same surface values as &, /3, y, and let e\ f\ . . . ? vj 
denote functions depending on them in the same way as #,/*, . . . , w depend on a, /3, y. 
Thus, by Taylor's theorem, 

»'-=^^-«)+f(/'-y)+|W-.)+S(»'-«)+>-»)+$(o'-o)+H, 

where H denotes the same homogeneous quadratic function of e } — <?, &c. that w is of e, 
&c. If for d—e, &c. we substitute their values by (111), this becomes 

, dw d(a! — a) , dw d(ot ! —ct) , dwd(tx!—oi) , , TT 

Multiplying by dxdydz^ integrating by parts, observing that oj— a, jS'—jS* y f -~y vanish 
at the bounding surface, and taking account of (109), we find simply 

^(w , —w)dwdydz=^Kdxdydz. (112) 

But H is essentially positive. Therefore every other interior condition than that speci- 
fied by (a, /3, y), provided only it has the same bounding surface, requires a greater 
amount of work than w to produce it : and the excess is equal to the work that would 
be required to produce, from a state of no displacement, such a displacement as super- 
imposed on (a, /3, y) would produce the other. And inasmuch as (a, /3, y) fulfil only the 
conditions of satisfying (110) and having the given surface values, it follows that no 
other than one solution can fulfil these conditions. 

68. But (as has been remarked by Professor Stokes to the author) when the surface 
data are of force, not of displacement, or when force acts from without, on the interior 
substance of the body, the solution is not in general unique, and there may be con- 
figurations of unstable equilibrium, even with infinitely small displacement. For in- 
stance, let part of the body be composed of a steel bar magnet ; and let a magnet be 
held outside in the same line, and with a pole of the same name in its end nearest to 
one end of the inner magnet. The equilibrium will be unstable, and there will be posi- 
tions of stable equilibrium with the inner bar slightly inclined to the line of the outer 
bar, unless the rigidity of the rest of the body exceed a certain limit. 

69. Eecurring to the general problem, in which the strains are not supposed infinitely 
small, we see that, if the solid is isotropic in every part, the function of A, B, C, #, 5, c 
which expresses w must be merely a function of the roots of the equation* 

(A~£ 2 )(B~^)(C-r)~<A~£ 2 )-^ . (113) 

which (that is the positive values of £) are the ratios of elongation along the principal 

* Thomson and Tait's < Natural Philosophy/ § 181 (11). 
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axes of the strain-ellipsoid. It is unnecessary here to enter on the analytical expression 
of this condition. For the case of A— 1, B— 1, C— 1, a, b, c, each infinitely small, it 
obviously requires that 

0> *)=(/> f)=(g, g) i (/» ff)=(g> *)=M) i 0> $)=& *)=(*> <0 ; ' 

(^, «)=(/, i)=(y,c?)=0; (J, (?)=((?, flj)=(a,i)=0; I . (114) 

fa *)=fa <>)=(/> <>)=(/> ^)=fa «)=fa *)=°- I 

Thus the twenty-one coefficients are reduced to three — 

(0, e), which we may denote by the single letter 3, 

\J ' 5/* " '* " ^ 5 ^' 

v^? <zj, „ ,, „ ,, n. 

It is clear that this is necessary and sufficient for ensuring cubic isotropy— that is to say, 
perfect equality of elastic properties with reference to the three rectangular directions 
OX, OY, OZ. But for spherical isotropy, or complete isotropy with reference to all 
directions through the substance, it is further necessary that 

$— 35=2^, .......... (115) 

as is easily proved analytically by turning two of the axes of coordinates in their own 
plane through 45° ; or geometrically by examining the nature of the strain represented 
by any one of the elements a, b, c (a "simple shear") and comparing it with the resultant 
of 6% and/= — e (which is also a simple shear). It is convenient now to put 

3-f2$=2m; so that %L=m+n 9 3&=m,-~n; ....... (116) 

and thus the expression for the potential energy per unit of volume becomes 

2 w = m (*+/+0)»+n^^ . . . (117) 

Using this in (108), and substituting for e, /, #, a, #, c their values by (111), we find 
immediately, for the equations of internal equilibrium, equations the same as (1) of 

§2. 

70. To find the mutual force exerted across any surface within the solid, as expressed 

by (26) of §• 22, we have clearly, by considering the works done respectively by P, Q, 
R, S, T, U (§ 21) on any infinitely small change of figure or dimensions in the solid, 

P=^, Q=fo E=*>, 8=% T= d ^i U=^. . . . (US) 

de df dg da do dc 

Hence, for an isotropic solid, (117) gives the expression (25) of § 21, which we have 

used above. 

71. To interpret the coefficients m and n in connexion with elementary ideas as to the 
elasticity of the solid; first let a=b=c=0, and e=f=g=^\ in other words, let the 
substance experience a uniform dilatation, in all directions, producing an expansion of 
volume from 1 to 1+& In this case (117) becomes 

S0=|(m— ^)S 2 ; 



m—^n 
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and we have 

Hence (m— ^n)l is the normal force per unit area of its surface required to keep any 
portion of the solid expanded to the amount specified by & Thus m~~^n measures the 
elastic force called out by, or the elastic resistance against, change of volume : and 
viewed as a coefficient of elasticity, it may be called the elasticity of volume. What is 

commonly called the " compressibility" is measured by 

And let next e=f=g=b=c=0 ; which gives 

w=j;na 2 ; and, by (118), S=na. 

This shows that the tangential force per unit area required to produce an infinitely 
small shear*, amounting to a, is na. Hence n measures the innate power with which 
the body resists change of shape, and returns to its original shape when force has been 
applied to change it ; that is to say, it measures the rigidity of the substance. 

[Note added, December 1863]. 

Since this paper was communicated to the Eoyal Society, the author has found that 
the solution of the most difficult of the problems dealt with in it, which is the determi- 
nation of the effect produced on a spherical shell by a prescribed application of force to 
its outer and inner surfaces, had previously been given by Lame in a paper published 
in Liouville's Journal for 1854, under the title "Memoire sur l'Equilibre l'elasticite 
des enveloppes spheriques. " In the same paper Lame shows how to take into account 
the effect of internal force, but does not solve the problem thus presented except for 
the simple cases of uniform gravity and of centrifugal force. The form in which the 
analysis has been applied in the present paper is very different from that chosen by 
Lame (who uses throughout polar coordinates) ; but the principles are essentially the 
same, being merely those of spherical harmonic analysis, applied to problems presenting 
peculiar and novel difficulties. 

* Thomson and Tait's ' Natural Philosophy,' § 171. 



